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Abstract 

We argue that the production of Mueller-Navelet jets at the LHC represents a convenient en- 
vironment to study gluon saturation and high-energy scattering in the presence of unitarity 
corrections. We show that, in a suitable range of transverse momenta for the produced jets, the 
cross-section for the partonic subprocess should exhibit geometric scaling. We point out that, 
in the presence of a running coupling, the cross-section for producing hard jets cannot be fully 
computed in perturbation theory, not even after taking into account the saturation effects: the 
non-perturbative physics affects the overall normalization of the cross-section, but not also its 
geometric scaling behavior. 
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1 Introduction 

Geometric scaling is one of the most important manifestations of the saturation physics 
in QCD at high energy [1-4], with striking consequences for the phenomenology. Identified 
first [1] in the HERA data, via a phenomenological analysis inspired by the idea of satura- 
tion [5] , geometric scaling has been soon after understood [2-4] as a property of the BFKL 
evolution [6] in the presence of saturation. For problems with a single (transverse) resolu- 
tion scale Q^, so like deep inelastic scattering (DIS) or inclusive single-particle production 
in hadron-hadron collisions, 'geometric scaling' means that, within a wide kinematical 
window, the cross-section scales as a function of the ratio Q"^ /Ql(Y). Here, QsiY) is the 
saturation momentum, i.e., the characteristic momentum scale for the onset of unitarity 
corrections in a collision in which the projectile and the target are separated by a rapidity 
gap Y. This scale grows rapidly with F, according to the BFKL evolution [2,3,7]. 

Most importantly, geometric scaling is not restricted to the saturation region at < 
Ql{Y), where the gluon occupation numbers are large and the scattering amplitudes are 
close to their unitarity limits, but it also extends over a relatively wide 'geometric scaling 
window' [2,3] at ^ Ql(Y), where the target is dilute and the scattering is weak, yet 
the scattering amplitudes 'feel' the effects of saturation, via the boundary condition at 
~ QliY). This scaling window, which with increasing Y is pushed towards larger and 
larger values of (because of the corresponding rise in Ql(Y)) and whose width is slowly 
increasing with Y (via the BFKL diffusion), is essentially the same as the validity range 
for the BFKL approximation [6] at high energy. Hence, the large-Q^ form of geometric 
scaling is a direct consequence of the BFKL dynamics precursory of saturation [2-4] and 
can be used to test the latter at the level of the phenomenology. 

So far, the most compeUing such tests have been performed on the HERA data at 
small values of Bjorken's x ^ Q'^/s, which are the data for which geometric scaling has 
been originally identified [1,8] (see also Refs. [9,10] for recent analyses, which include the 
diffractive data). Namely, one found that, within the whole small-x domain at HERA, 
i.e., for X < 0.01 and Q'^ < 450 GeV^, the DIS cross-section can be well approximated 
by a scaling function: a{x, Q^) ^ (t{Q'^/QI{Y)) with Ql{Y) ~ e^'^ and Y = ln(l/x). (For 
comparison, the proton saturation momentum at HERA is estimated in the ballpark of 1 
GeV.) Remarkably, the value A ~ 0.3 for the 'saturation exponent' emerging from these 
analyses is in rough agreement with its perturbative calculation [7] from the next-to- 
leading order BFKL equation [11,12]. More detailed analyses [13-16], combining BFKL 
dynamics and unitarity corrections within the framework of the dipole picture, have shown 
that the HERA data are consistent with some of the hallmarks of the BFKL evolution, like 
its characteristic 'anomalous dimension', or the violation of geometric scaling via BFKL 
diffusion. Similar parametrizations for the dipole cross-section, with the parameters fixed 
through fits to the HERA data, have been used [17-19] to describe particle production 
in deuteron-nucleus collisions at RHIC, with some success in explaining the 'high-pr 
suppression' in the nuclear modification factor at forward rapidities. 

However, given the kinematical limitations inherent in the experiments at HERA and 
RHIC, the previous phenomenological studies of BFKL physics, geometric scaling, and 
saturation cannot be viewed as definitive. The situation should be more favorable in 
this respect at LHC, where the higher available energies and the experimental setup 
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should offer larger rapidity gaps to the BFKL evolution. For instance, in forward particle 
production at LHC one could measure values of x as small as a; ~ 10~^ in the 'target' 
proton wavefunction for a produced jet with transverse momentum k_[_ ~ 10 GeV. With 
this kinematics, the jet should explore the geometric scaling window of the target proton, 
with interesting consequences, e.g., for the nuclear modification factor [20,21]. 

Another interesting process that was proposed to test the BFKL dynamics and which 
could be measured at LHC under favorable conditions is the production of Mueller-Navelet 
jets [22] (see also Refs. [23-31] for various theoretical studies and [32,33] for experimental 
searches at the Tevatron). This is a pair of jets separated by a large rapidity gap Y which 
should favor the BFKL evolution of the cross-section for the partonic subprocess. For 
sufficiently large values of F, saturation effects (in the form of unitarity corrections to 
the partonic scattering) should become important, as already emphasized by Mueller and 
Navelet in their original proposal [22]. However, with the exception of a few, preliminary, 
phenomenological studies [34,35], such effects have been left out in previous studies of 
the Mueller-Navelet jets, which focused on the energies at the Tevatron. In particular, 
the modification of the BFKL dynamics by saturation and the phenomenon of geometric 
scaling have never been addressed in this context. These are the aspects that we would 
like to focus on in what follows. 

Our first observation is that the Mueller-Navelet process is particularly favorable to 
study saturation physics. Unlike in DIS, where the gluon evolution inside the proton starts 
at the 'soft' scale Aqcd ~ 250 MeV and thus requires a relatively large rapidity evolution 
before it develops a hard saturation momentum, in the context of Mueller-Navelet jets 
this evolution starts with the 'hard' scale set by the transverse momentum > 10 GeV 
of one of the two jets. Hence, the subsequent evolution with Y should rapidly produce a 
system with very high gluon density around the position of the jet, i.e., with a very large 
local saturation momentum. The counterpart of that is that the dense region occupies 
only a small area ~ 1/^1 in impact parameter space, and looks like a 'dense spot'. 

To be more specific, recall that one needs a rapidity evolution Yq ~ (l/tup) ln(l/a^), 
with ujf the BFKL intercept, before a small hadronic system, so like a dipole or a high- 
momentum parton, reaches saturation on the resolution scale set by its own size, or 
transverse momentum [36]. A leading-order estimate for ujv would yield Iq — 5 (for 
as = 0.2, as appropriate for a 10 GeV jet), but this is probably too optimistic. A more 
realistic estimate, using the NLO BFKL intercept [12], is Yq ~ 8, which is fully within the 
reach of LHC. This means that, when producing a pair of Mueller-Navelet jets separated 
by a rapidity gap y = 8 at the LHC, then one of the jets will 'see' the other one as a high- 
density gluonic system (a 'color glass condensate') with a saturation momentum Qs{Y) ~ 
10 GeV. Moreover, every additional unit of rapidity will make this saturation scale even 
harder, according to Ql(Y) ~ Qlexp[X(Y — Yq)]. Such values for Qg are considerably 
higher than those that could ever be achieved in a proton, or nuclear, wavefunction at 
LHC energies, even for the most forward collisions. 

Furthermore, the evolution towards saturation should favor the transverse momentum 
dissymmetry between the two jet^. The typical transverse momentum of a gluon within 



^ We recall that in lowest-order perturbation theory the two jets come out with equal and oppo- 
site transverse momenta, because of momentum conservation. Hence, any momentum asymmetry 
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the wavefunction of an evolved jet is the respective saturation momentum Qs{Y), which 
for large Y is as hard as, or even harder than, the original parton which comes out as 
a jet. Therefore, the other jet can easily be produced with a very different transverse 
momentum, because the momentum imbalance can be compensated by 'inclusive' gluons 
from the 'evolved' jet wavefunction. In fact, we expect a larger momentum asymmetry 
in the presence of saturation than from pure BFKL evolution: unlike the latter, which 
proceeds symmetrically towards soft and hard momenta, the evolution in the presence of 
saturation is biased towards large transverse momenta (larger than the saturation scale). 

The most interesting kinematical situation for our subsequent analysis is precisely when 
the two jets are well separated in transverse momentum, say, ki± ^ k2±- More precisely, 
we shall be interested in configurations where the harder jet has a transverse momentum 
comparable to, or even larger than, the saturation momentum that would be generated 
by the evolution of the softer jet over the rapidity gap Y: ki± > Qs{Y), where Qs{Y) 
is implicitly a function of k2j_- Under these circumstances, we shall see that the partonic 
cross-section exhibits geometric scaling within a wide kinematical window. That is, for 
given /c2±, the cross-section scales as a function of the ratio kfj_/Q'^{Y). On the other 
hand, the cross-section is small, of order l/k2j_, because of the small size of the dense 
spot, as alluded to above. While our conclusions may look natural, given the kinematics 
and the similarity with other problems like DIS, our analysis appears to reserve some 
difficulties and surprises. 

The first difficulty refers to the inclusion of saturation effects and unitarity corrections 
in the cross-section for Mucllcr-Navclct jets. This in turn requires two steps: (i) a fac- 
torization formula which is general enough to allow for unitarity corrections, and (ii) the 
calculation of the ingredients which enter this factorization formula within the framework 
of high-density QCD (i.e., from the solutions to the non-linear evolution equations which 
generalize the BFKL equation to the region of high gluon density). 

Concerning step (i). wc shall proceed in a heuristic way, by generalizing, in Sect. 2, a 
known formula for single-jet production in the presence of unitarity corrections [37-41]. 
This leads us to a generalization of the standard /cr-factorization for Mueller-Navelet 
jets, which has been already presented in Ref. [34,35,41], and in which the BFKL Green's 
function is replaced by the total cross-section for the scattering between two effective 
color dipoles. Unlike the quark-ant iquark dipole familiar in the context of DIS (see, e.g., 
[5, 13-16]), which is a physical fluctuation of the virtual photon, the dipoles that enter 
our factorization for Mueller-Navelet jets are merely mathematical constructions, which 
appear in the calculation of the cross-section and are built with one parton in the direct 
amplitude and another parton in the complex conjugate amplitude. The dipole-dipole 
cross-section is written in coordinate space, as appropriate for the inclusion of unitarity 
corrections (multiple scattering) in the eikonal approximation. The transverse momenta 
ki± and k2± of the produced jets are then fixed via a double Fourier transform from the 
dipole sizes. 

The second step, i.e., the calculation of the dipole-dipole cross-section within high- 
density QCD, turns out to be particularly subtle. Since the scattering involves two systems 



between the jets is a signal of a multiparticle final state, as produced in particular by the BFKL 
evolution. 
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(dipoles) which start by being dilute at low energy, it seems that we cannot rely on the 
standard Balitsky-JIMWLK, or BK, equations [42-45], which apply only to dense-dilute 
scattering. Instead, one should use the more general, 'Pomeron loop', equations [46-48], 
which also allow for particle number fluctuations in the course of the evolution. From 
the correspondence with statistical physics [49], and also from the numerical simulations 
of simple models inspired by QCD [50, 51], we know that, with a fixed coupling, the 
effects of the fluctuations are truly crucial: with increasing energy, they rapidly wash 
out both the BFKL approximation and the 'geometric scaling' behavior predicted by the 
BK equation [2-4]. However, a very recent numerical analysis [52] has shown that the 
fluctuations are strongly suppressed by the running of the coupling, in such a way that 
their effects remain negligible for all energies of practical interest. 

With the philosophy that the running-coupling case is the only one of fundamental 
interest for real QCD, in what follows we shall perform a 'mean field' type of analysis, 
based on BK equation, for both fixed and running coupling. The fixed-coupling analysis, 
as developed in Sects. 3 and 4, turns out to be rather straightforward: The dipole-dipole 
cross-section, as obtained from approximate solutions to BK equation [2-4], exhibits 
geometric scaling for suitably chosen dipole sizes. After a Fourier transform, this scaling 
property gets transmitted to the partonic core of the Mueller-Navelet cross-section, for 
appropriate transverse momenta of the two jets. 

The running-coupling case, that we shall treat in Sect. 5, is both more interesting and 
more subtle. First, it might look inconsistent to include the running of the coupling, but 
at the same time ignore other next-to-leading order corrections in perturbative QCD. 
But it turns that the running of the coupling plays a special role in the context of the 
high-energy dynamics: because of it, all the other perturbative corrections die away in 
the high-energy limit [7, 53] . Indeed, the evolution towards saturation is controlled by 
momenta k± ~ Qs(Y); so, the relevant value of the coupling is asiQ^Y)), which decreases 
with Y, and therefore so do the perturbative correction^, whose strength is proportional 
to as{Ql{Y)). Moreover, the running of has qualitative consequences which modify the 
high-energy evolution in depth. We have already mentioned its role in suppressing the 
particle-number fluctuations. This is related to a more general property of the running 
of the coupling, which is to slow down the evolution towards saturation [2,3,7]. Another 
manifestation of this property is in the growth of the saturation momentum with Y: for 
sufficiently large Y, and with a running coupling, \nQl grows like y/Y, and not like Y. 

An important consequence of the running of the coupling, which is not specific to the 
high-energy problem, but has dramatic consequences for it, is the fact that it introduces 
an intrinsic scale in the problem — the 'soft' scale Aqcd — , thus breaking down the 
conformal invariance of the leading-order formalism. With increasing energy, this scale 
progressively replaces within the saturation momentum any other scale introduced by 
the initial conditions at low energy, so like the target dipole size. Accordingly, for suffi- 
ciently high energy, the saturation momentum becomes independent of the target size [54] . 
This has important consequences for the high-energy evolution in general (e.g., it implies 
that a large nucleus is not more dense than a proton at very high energies), and for the 



^ Of course, these corrections may be numerically important for the phenomenology at LHC, 
but here we focus on the dominant asymptotic behavior, for simplicity. 
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Mueller-Navelet process in particular: it implies that the perturbative calculation of the 
Mueller-Navelet cross-section breaks down, even if the jet transverse momenta are re- 
stricted to be hard. The precise argument in that sense will be developed in Sect. 4, but 
here we would like to emphasize that this argument is in fact very general (and hence also 
very robust): it reflects the fact that, with running coupling, the dipole-dipole scattering 
amplitude at a fixed impact parameter is independent of the target dipole size R (rather 
than dying away as an inverse power of R, as it would happen in the fixed-coupling for- 
malism, by conformal invariance). Accordingly, the dipole-dipole cross-section, which is 
obtained by integrating the amplitude over all impact parameters, is proportional to R^, 
and hence it strongly favors large dipole fluctuations. Without the non-perturbative cut- 
off introduced by confinement, the partonic cross-section would be controlled by dipole 
fluctuations with arbitrarily large size. As we shall argue in Sect. 4, the ad-hoc intro- 
duction of a non-perturbative cutoff on the dipole sizes affects the normalization of the 
total cross-section, but not also its property of geometric scaling (which merely refers to 
the functional dependencies of the cross-section upon the rapidity gap Y and upon the 
transverse momentum ki_i_ of the hardest jet). 



2 Forward jets with unitarity corrections 

Although our main interest here is in the production of a pair of (Mueller-Navelet) jets, 
it is instructive to start our presentation with the case of a single jet, for which the high- 
energy factorization in the presence of unitarity corrections is more firmly established. 
This will also allow us to introduce the physics and the theoretical description of the 
unitarity corrections in a simpler setting. The two-jet problem will then be easier to 
explain, by analogy. 

For definiteness, we focus on jets initiated by gluons (quarks will be added later on), 
and thus consider the cross-section for inclusive gluon production at forward rapidity in a 
hadron-hadron collision at high energy. By 'forward rapidity' we mean that the produced 
gluon carries a sizeable fraction x ~ 0(1) of the longitudinal momentum of one of the 
incoming hadrons (the 'projectile'), so that there is a large rapidity gap Y = Yq — y 
between this produced gluon and the other hadron (the 'target'). Here, Yq = \n{s / M1M2) , 
with s the invariant energy squared and Mi^2 the masses of the participating hadrons, is 
the rapidity gap between the projectile and the target, and y = ln(l/x) +ln(A;_L/Mi) is the 
(relatively small) rapidity separation between the produced gluon, which has transverse 
momentum k±, and the projectile. Alternatively, one could trade the rapidity gap y for 
the pseudo-rapidity rj of the produced jet in the laboratory frame; e.g., if the lab frame 
coincides with the hadron center-of-mass frame, so like at LHC, then x = {k±/ \/s) e"^ , 
with ?7 > for forward jets (see also Fig. [Ha). 

Under these circumstances, and in the leading-order formalism of perturbative QCD 
(meaning, in particular, that the coupling is fixed), the cross-section for gluon production 
can be expressed in a '/c^-factorized' form, which is formally similar to, but more general 
than, the corresponding factorization used in the context of the BFKL physics [6]. Namely, 
the /cr^factorization is now extended towards the high-energy regime where unitarity cor- 
rections (multiple scattering, gluon saturation) become important. The general respective 
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formula can be found in the literature (see, e.g., Refs. [37-41,55]). Here, we shall need 
only a simpler form of it, valid when the transverse momentum k_i_ of the produced gluon 
is large enough — much larger than the typical momentum transferred from the projectile 
to this gluon (see below for a more precise condition). We then have 

] PT^JX 1 ^ 

= Si5-G.(.,*i)/dV e— V!a,„,r(r,Y), (2.1) 

where xGp{x, k\) is the gluon distribution in the projectile (P) on the resolution scale of 
the jet (i.e., the number of gluons with longitudinal momentum fraction x equal to that 
of the jet, and with transverse momenta p\ < k'jj. Furthermore, CF(^gg)T{r, Y) is the total 
cross-section for the scattering between a gluonic dipole (a gg pair in a color singlet state) 
with transverse size r and the hadronic target (T), for a rapidity separation Y . The gg 
dipole here is the effective dipole made with the produced gluon in the direct amplitude 
(located at transverse coordinate x) and the corresponding gluon in the complex conjugate 
amplitude (located at y). The gluon transverse momentum fc_L = |fc| in the final state is 
then fixed via the Fourier transform from r = x — y to k. 

In the single-scattering approximation to the dipole-target cross-section, the Fourier 
transform in Eq. (12. ip yields the usual 'unintegrated' gluon distribution in the target 
wavefunction (evaluated in the BFKL approximation) times a constant of order a^. We 
then recover from Eq. (12.11) the traditional /cT^factorization. But Eq. (12. ip remains valid 
also very large values ofY, where the unitarity corrections to the dipole scattering become 
important and the BFKL approximation ceases to apply. 

The physical interpretation of the unitarity corrections is most transparent in the 
'target infinite momentum frame', where the dipole has relatively low energy while the 
target carries most of the total rapidity Y. Then the target wavefunction has evolved into 
a 'color glass condensate' (CGC) — a system with high gluon density characterized by 
a hard intrinsic scale, the saturation momentum Qs(Y), which grows rapidly with Y and 
separates between 

• a high density region at low momenta p± < Qs{Y), where the gluon occupation numbers 
are large, ~ 0{l/as), but 'saturated' (they do not rise with the energy anymore), and 

• a low density region at high momenta p± > Qs{Y), where the occupation numbers are 
low, but rapidly growing with Y, via the BFKL evolution. 

For sufficiently high energy and/or large dipole sizes, such that r > l/Qs{Y), the dipole 
will undergo multiple scattering off the CGC. The 'unitarity corrections' refer to this 
multiple scattering, as well as to the saturation effects in the target gluon distribution, i.e., 
to all the non-/mear effects which reduce the gluon density and enforce the unitarity bound 
on the scattering process. Such effects are resummed — within the eikonal approximation, 
and within the limits of the LO formalism — by the Balitsky-JIMWLK equations, which 
in particular determine the dipole scattering amplitude in this high energy regime ^ I. The 
corresponding solution will be further described in the next sections. Here, it suffices to 

^ As mentioned in the Introduction, these equations neglect the particle number fluctuations, 
which would be important in the context of the fixed-coupHng evolution [46,49], but which are 
suppressed by the running of the coupling [52]. 
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(a) 



(b) 



Fig. 1. The inclusive production of a single, forward, jet (left) and of a pair of Mueller-Navelet 
jets (right ) in the presence of unitarity corrections. 

notice that the factorization (12.11) is correct when Y is so large that Ql{Y) ^ ^qcd (i^ 
order for the perturbative approach to unitarity corrections to be vahd), and for relatively 
'hard' jets, whose transverse momenta are not much smaller than Qs{Y) — the precise 
condition being kj_ 3> Qs(Y) Aqcb- Indeed, under these circumstances, the contribution 
shown in Eq. (12.11) is enhanced with respect to the other, missing, contributions by the 
large logarithm \n(^kj_/ Aqqj)) (via the gluon distribution xGp{x, k'j_)). 

We are now prepared to present the corresponding formulae for the Mueller-Navelet 
jets. This is a pair of jets produced in a high-energy hadron-hadron collision such that 
each jet carries a relatively large fraction Xj, i = 1, 2 of the longitudinal momentum of its 
parent hadron. Accordingly, each jet is relatively close in rapidity to its respective parent, 
so there is a large rapidity gap Y ^ 1 between the jets: Y = \rL(xiX2S / ki±k2±) , where 
ki± and k2± are the jet transverse momenta (see Fig. [Hb). When the scattering is viewed 
in the hadron center-of-mass frame, the momentum fractions Xi and the rapidity gap are 
determined by the pseudo-rapidities rji of the two jets, according to 



The typical kinematics for Mueller-Navelet jets is such that rji is large and positive, while 
ri2 is large and negative. 

The differential cross-section for Mueller-Navelet jets has been rigorously computed 
[22] (within the LO formalism, once again) only at the level of the BFKL approximation, 
which ignores unitarity corrections. The corresponding result is the expected generaliza- 




k2± . 
X2 — —p= e 



Y 



(2.2) 
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tion of the corresponding single-jet cross-section — the BFKL version of Eq. (12.11) — 
which is symmetric w.r.t. the two jets. In view of this, and of the symmetry of the prob- 
lem, it has been conjectured [34, 35, 41] that, after including the unitarity corrections, 
the cross-section for Mueller-Navelet jets should be given by the properly symmetrized 
version of Eq. (12.11) , that is 



X 



J dVie-"=-Y d'rse-^-^^ V,\V^^a(,,)(,,)(ri,r2,F), (2.3) 



where, for definiteness, we have chosen the incoming hadrons to be protons, so like at 
LHC. o'{gg){gg){'fi, ^2, ^) IS the total cross-scctiou for the scattering between two (effective) 
gluonic dipoles with the indicated transverse sizes and separated by a rapidity gap Y. 

The cross-section (12. 3p is generally a function of ki±, k2±, and the relative angle 
between the vectors ki and ^2 (the azimuthal angle between the two jets). Although 
interesting in view of the phenomenology (and largely studied in the context of the BFKL 
approximation [23-31]), the azimuthal correlations represent a subleading effect at high 
energies, and thus are irrelevant for our present study of the unitarity corrections. So, in 
what follows we shall average over 0, which is tantamount to replacing the dipole cross- 
section in (12. 3p by (7(^gg)(^gg){ri,r2, Y) — the corresponding cross-section averaged over the 
relative orientations of the two dipoles. 

It is furthermore convenient to introduce some hard momentum cutoffs Qi and Q2 
(corresponding to the experimental fcy-cuts) and compute the cross-section for producing 
two jets with transverse momenta ki± > Qi and k2i_ > Q2 and with given longitudinal 
momentum fractions Xi and X2- That is, 

(These integrations automatically implement the average over 0.) A priori, the inte- 
grations are complicated by the /c_l -dependencies of the gluon distributions in (12.31) . 
and by the one implicit in the dipole-dipole cross-section, via the rapidity gap Y = 
ln{xiX2s/ki±k2±). Note however that for given xi, X2 and s, the maximal rapidity gap 
^max = Iti{xiX2s/QiQ2) is attained for the threshold momenta ki± = Qi and k2± = Q2- 
Hence the dominant contribution, in the sense of the leading-logarithmic approximation, 
is obtained by replacing Y — ^ l^max within o-(gg)(^gg). Moreover, the integrand in (12. 4p is 
rapidly decreasing at very large values for ki± and k2± — this can be checked, e.g., by us- 
ing the BFKL approximation for the dipole-dipole cross-section [22] — , which enables us 
to replace G(xj, kf^) G(xj, Qf) in the slowly varying gluon distributions. The remaining 
integrations can be easily performed, with the final result 

^i pp^jxj 1 

G{xi,Q\)G{x2,Ql) 



00 00 



da;ida;2 647r^ 

X J dri J dr2QiJi{Qiri)Q2Ji{Q2r2)cr{gg)(gg)iri,r2,Y). (2.5) 
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Here and from now on it is understood that Y = \n{xiX2s/QiQ2)- So far, we have con- 
sidered only gluon jets, but quarks or antiquarks jets can be similarly included: when the 
jet i, with i = 1,2, is initiated by a quark with flavor /, we have a formula similar to 
Eq. (12. 5p in which the gluon distribution G{xi, Qf) is replaced by the quark distribution 
g/(xj, Qf) (or g/(xj, Qf) for an antiquark), and the corresponding dipole within o-(^gg)(^gg) is 
replaced by a dipole made with a quark and an antiquark. We thus encounter three types 
of dipole-dipole processes: {gg){gg), {,gg){(l'i), and {qq){qq), whose cross-sections differ at 
most through color factors (see below). 

As a consistency check of our above factorization of the Mueller-Navelet cross-section, 
cf. Eq. (12.31) or (12.51) . let us now verify its BFKL limit. For two gluonic dipoles, the BFKL 
cross-section (averaged over angle) read£3 



BFKL 



Cf J 2m 7^(1-7) 



n 



X{1)Y 



(2.6) 



where we use the standard representation for the BFKL solution in Mellin space (see, 
e.g., [2,3]). This result is symmetric under the exchange ri ^ r2 of the two dipoles, as 
it can be checked by using the property xil) = ~ 1) of the BFKL characteristic 
function. When one or both of the gluonic dipoles are replaced by fermionic (qq) ones, the 
expression in Eq. (12. 6p must be multiplied by a factor Cp/Nc for each such a replacement. 
Substituting (12.61) into (12. Sp . making use of the Bessel function integration formula 

^/ n + l + /3 \ 

dxx^Ux) = 2^ )^^l_J , (2.7) 



^' 2 

and summing up over all types of partons, we obtain (with a = agNc/ vr) 

da^^-'''' p MN^al rd7(Q|/Q|r , 
= —T- 752 / 77- —r. -r exp \ax{r()Y\, {2X 



dxida;2 



where Fgfj involves contributions from quarks, antiquarks, and gluons, with appropriate 
color factors (the sum over the quark flavors is kept implicit): 

-^cff = „. . /cff (3^1; Ql) /eff(^2; Q2) 

647r* 

/eff (x, Q2) = G{x, Q') + ^ [qix, Q') + qix, Q')] . (2.9) 

As anticipated, Eq. (12. 8p is in precise agreement (including the normalization) with the 
corresponding result in Ref. [22]. 



^ The appearance of the Casimir Cf for the fundamental representation in a cross-section 
pertinent to gluons alone may look surprising. In reality, this has been generated via the identity 
N^/{N^ — 1) = Nc/2Cf, where all the factors arise from the gluon color algebra. 
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3 Fixed coupling case: the dipolar cross— section 



In addition to taming tlie BFKL growtli of tlie dipole scattering amplitude, in compli- 
ance with the unitarity bound, the non-linear effects encoded in the Balitsky-JIMWLK 
(or BK) equations have also an interesting consequence for the functional form of the 
amplitude in the transition region from weak to strong scattering: within a rather wide 
kinematical region, whose width is increasing with Y, this amplitude shows geometric 
scaling [2,3], i.e., it depends upon the size r of the projectile dipole and upon Y only 
via the dimensionless variable r = r'^Q'j.{Y), with Qs{Y) the saturation momentum of 
the target evolved up to rapidity Y. It is then tempting to conjecture that this scal- 
ing property should transmit from the dipole-dipole amplitude to the cross-section for 
Mueller-Navelet jets, via the convolutions in Eq. (12. 5p (within a suitabe range of values 
for the variables Qi, Q2, and Y). This is a simple argument in that sense: due to the 
presence of the rapidly oscillating Bessel functions in the integrand of Eq. (12. 5p . one ex- 
pects the integrals there to be dominated by values ri ~ l/Qi and r2 ~ 1/Q2- If this is 
true, then one can choose Qi and Q2 (for a given Y) in such a way that the dipole-dipole 
cross-section, and hence the dijet cross-section, are in the geometric scaling window. As 
we shall later discover, via explicit calculations, this simple argument is indeed correct in 
the case of a fixed coupling, but not also for a running coupling. 

We start with the fixed-coupling case, i.e., the LO formalism. Within the eikonal 
approximation, the dipole-dipole cross-section is computed as (for dipoles made with 
partons in a generic color representation) 



where Tdd(T'i, ^2, b, Y) is the scattering amplitude for two dipoles with transverse sizes ri 
and r2, relative impact parameter b, and rapidity separation Y. (The average over the 
relative angle between Vi and r2 is implicit here and from now on.) We use conventions in 
which the S'-matrix is written as = 1 — T, where T is taken to be real, as appropriate 
for the dominant behavior at high energy. In general, the dipole-dipole amplitude and 
cross-section are, of course, symmetric under the exchange of the two dipoles, but their 
approximate forms that we shall derive below are valid only when one of the dipoles 
is much smaller than the other one. It is then convenient to introduce the notations 
r = min(ri, r2) and R = max(ri, r2), and refer to the small (large) dipole as the 'projectile' 
(respectively, the 'target'). Also, as explained in Sect. [2], it is convenient to visualize the 
evolution with increasing energy as gluon evolution in the 'target' (the larger dipole). 

The unitarity of the S'-matrix implies T < 1, with the upper bound T = 1 (the 'black 
disk limit') describing a situation where the scattering occurs with probability one. This 
constraint is indeed obeyed by the solution T to the Balitsky-JIMWLK (or BK) equations, 
which, moreover, appears to saturate the black disk limit T = 1 for sufficiently large Y. 
But the cross-section (13. ip can rise indefinitely with Y, even after the 'black disk' limit 
has been reached at central impact parameters, because the gluon distribution in the 
target keeps expanding towards larger impact parameters, due to the non-locality of the 
BFKL evolution. This radial expansion of the gluon distribution is however much slower 
than its evolution towards the black disk limit at a fixed value of b. In particular, we do 




(3.1) 
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not expect this expansion to be essential at the LHC energies. Therefore, a target which 
at y = starts as a single dipole of size R, evolves with Y towards blackness, first, on its 
own scale R, then, on smaller and smaller scales, without significantly expanding towards 
larger size^U. Accordingly, when this evolved target is probed by a projectile dipole with 
size r <^ R, the amplitude Tdd{r, R, b, Y) is negligibly small when the two dipoles have no 
overlap with each other {b ^ R). The typical impact parameters which contribute to the 
cross-section are such that b <^ R, and for them the amplitude is roughly independent of 
b. These considerations motivate the following approximation to the dipole-dipole cross- 
section (13.11) . valid when r <ti R 

(Tdd (r, R, Y) ~ 2rf Tdd(r, R, Y) , (3.2) 

where Tdd{r, R, Y) is independent of b and satisfies the unitarity bound Tdd < 1- For this 
amplitude, we shall use approximate solutions to the BK equation with fixed couphng. 
Specifically, the saturation momentum Qs{R,Y) is defined by the condition 

Tdd{r,R,Y) = K for r = l/g,(i?, F) , (3.3) 

where k < 1 is a number of order one (its precise value is irrelevant to the accuracy of 
interest). For r > l/Qs we have Tdd ~ whereas for r <^ 1/Qs the amplitude is 

small, Tdd 1, and approximately given by the following, universal, function 

T«(.«.yO- (ln-i^)(.^Ojpexp{-!^M} (3.4) 

(up to a normalization factor), with the saturation momentum 

QliR,Y) = (Ca,2)^/^^-^e^-^. (3.5) 

The various exponents which appear in these formulae are pure numbers determined by 
the BFKL characteristic function (see Refs. [2,3] for details). Specifically, 7^ ~ 0.63 
(1 — 7s ~ 0.37 is the BFKL anomalous dimension at saturation), Dg ^ 48.5 plays the 
role of a diffusion coefficient, and Xg ~ 4.88 is the saturation exponent. Note that rapid 
growth of the saturation momentum with Y. Its dependence upon the target size R could 
have been anticipated from dimensional arguments. 

Eq. (13.41) is universal in the sense that the dependence upon the initial conditions at 
low energy is fully encoded in the value of the saturation momentum. But the latter is, 
of course, process-dependent: it depends upon the size R of the target, and also upon 
the color representations of the partons making up the two dipoles participating in the 
collision (via the coefficient C which here is left unspecified). For instance, the value of 

^ This property is not correctly encoded in the Balitsky-JIMWLK, or BK, equations, which 
rather predict a rapid radial expansion of the black disk, because of the long-range tails (in 
b) of the perturbative gluon distribution [56]. However, in real QCD we expect this tails to be 
cut-off by confinement, with the effect that the radial expansion is drastically slowed down, in 
compliance with Froissart bound [57]. 
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C corresponding to the process {gg){qq) is larger by a factor Nc/Cp than that for the 
process {qq){qq). 

The prefactor involving in the expression fl3.5p for the saturation momentum reflects 
the fact that the scattering amplitude starts at order in perturbation theory, hence 
one needs some non-trivial rapidity evolution Yq before the dipole becomes 'black' on 
the resolution scale fixed by its own size: T^ir = R,R,Yo) = k. This condition yieldJ^ 
Yq ~ (l/7<jAsa) ln(K/a^), and then Eq. fl3.5p can be rewritten in such a way to exhibit the 
rapidity excess beyond Yq: 

\sa{Y-YQ) 

QliR,y) = • (3.6) 

Eq. (13. 4p is valid when oY ^ 1 and in a rather wide range of values for r, namely for 
1 <C \a{l/r'^Ql) < caY , with c ~ 0{1). In particular, within the more restricted window 

1 



1 < ^^-^ < ^^DsaY, (3.7) 



the last, Gaussian, factor in Eq. (13.40 can be ignored, and then the amplitude shows 
geometric scaling, as anticipated. Note that, when oY ^ 1, this scaling window is quite 
wide, especially since the parameter Dg is numerically large. 



4 Fixed coupling case: Mueller— Navelet jets 

We now have all the ingredients to compute the cross-section for Mueller-Navelet jets in 
the presence of unitarity corrections and for fixed coupling. We shall focus on the range of 
values for the momentum cutoffs Qi and Q2 in which we expect geometric scaling. Namely 
we choose Qi ^ Q2 in such a way that Q\ > Q\ e'^^, with the compact notations A = A^o; 
and y = Y — Yq. That is, y is the rapidity excess introduced in Eq. (13. 6p . (Throughout 
this paper, we assume Y > Yq.). Therefore, Qi is larger, but not much larger, than the 
saturation scale that would generated by a target dipole with size ~ 1/(52 after a rapidity 
evolution Y . As we shall shortly check, under this condition the convolutions in Eq. (12. 5p 
are indeed dominated by values for ri and r2 within the scaling window (13. 7p . 

To simplify the calculation, we therefore keep only the scaling piece in the dipole-dipole 
amplitude (13.41) in the weak scattering regime. We therefore replace Eqs. (I3.2p - (l3.4p by 
the following, piecewise, approximation to the dipole-dipole cross-section (we recall the 
notations r = min(ri, r2) and R = max(ri, r2)) 

(4.1) 



* In the Introduction, this was written as Yq ~ (l/wp) ln(l/a^), which is essentially the same, 
since 7sAsa plays the role of the intercept at the saturation saddle point, as manifest in Eq. ()3.4p . 



7s 



for 
for 
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where we have ignored the slowly varying logarithm in Eq. fl3.4p . The precise normalization 
of the cross-section is not an issue here, as we are merely interested in its functional 
dependencies. 

Because of the symmetry of the above cross-section under ri ^ r2, when evaluating the 
double integral in Eq. (12.51) it is enough to consider the case ri < r2. Then the contribution 
from the other region ri > r2 can be simply obtained by letting Qi ^ Q2 in the result of 
the first cas^. Thus, substitution of (14.11) into (12. 5p gives 

— = Fe-'^^y / driQiJi(Qiri)r^= / drs Q2 Ji(Q2r2)r2'-'^» 

dxida;2 J J 



+ F JdnQMQ.n) J dr2Q2MQ2r2)rl + {Qi ^ Qs}, (4.2) 

ri 

where for the time being we consider the gluon jets alone (hence, the overall factor F 
includes the gluon distributions, together with other numerical factors); the quark jets 
will be added later on. Note that the exchange Qi Q2 should not be done inside F, 
but only in the result of the integration. Clearly, the two explicit terms in the r.h.s. of 
the above equation arise from the corresponding pieces in (14.11) . It is convenient to change 
variables by letting Ui = QiVi and U2 = Q2^2- Then the above equation becomes 

^^pp^JXJ ^ /Q2^\yys 



^ ^ aui 

^ cxD aui 

+ ^Q2j d^l Ji(Mi) J dU2ulJi{u2) + {Ql^Q2}, (4.3) 
bui 

where we set a = {Q2/Q1) exp(A?//2) and b = Q2/Qi- In what follows we shall show that 
the only non-zero term in the r.h.s. of Eq. (14.31) is the first one, and moreover this term 
shows geometric scaling. 

The fact that the second term in (14. 3 p vanishes could have been anticipated, since this 
term arises from the saturation piece in Eq. (14. ip . which in turn depends only on one 
of the two variables ri and r2. Hence, if we return to the unintegrated version of the 
cross-section, Eq. (12.31) . it becomes obvious that this saturation piece goes away by the 
successive action of the two Laplacians V^^V^^- I^ ^ little bit more difficult to see the 
corresponding cancelation in (14.30 . Performing the integration over U2 we obtain for this 
second term under consideration 

1 °° 

F-^ I duiulJiiui) [0^2(0^1) - 6V2(&Mi)l . (4.4) 



To evaluate this integral, let us choose n = 1 m the completeness formula 



^ Given our choice that Qi ^ Q2-, one could anticipate that the dominant contribution comes 
from the region ri <C 7^2. For completeness, we shall nevertheless consider the region ri > r2 too. 
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duuJn{u)Jn{au) = 5{a — 1), (4.5) 



then differentiate this identity with respect to a and use the fact that J[{x) = Ji{x)/x — 
J2(x). We thus obtain 



duu^Ji{u)J2{au) = 6{a - 1) - 6' {a - 1), (4.6) 



which provides the result for the integrals in Eq. (14. 4p . Namely, since we are interested in 
momenta such that Qi > Q2 6xp(A?/), we see that a < 1 and 6 < 1, hence both terms in 
Eq. (14. 4p vanish, as anticipated. 

Let us now turn to the calculation of the first term in (14. 3p . Putting aside the prefactor, 
the remaining double integration, let us call it h{a), can be simplified by differentiating 
with respect to a. We have 

oo 

h'{a) = -a2~27. J duu^Ji{u)Ji{au). (4.7) 



The above integral can be again recognized as the derivative of a known integral: by 
differentiating (14.61) with respect to a and using J2{x) = Ji(x) — 2J2{x)/x, we find 

oo 

J duu^Ji{u)Ji{au) = -S'{a - 1) - S"{a - 1). (4.8) 



It is now straightforward to obtain h{a) by integrating over a. Integrating by parts to get 
rid of the derivatives acting on the 5-function and using h{oo) = 0, we obtain 

/i(a) = 47,(l-7,)e(l-a). (4.9) 

Recalling that a < 1 in the kinematical region of interest, we see that the step function is 
equal to 1. That is, the first (double) integral in the r.h.s. of (14. 3 p is independent of a so 
long as a < 1. 

By a similar argument, it is now easy to see that the term obtained by exchanging 
Qi ^ Q2 (that is, the term coming from the region ri > r2) is equal to zero: indeed, the 
corresponding contribution would be proportional to 0(1 — 1/a). 

Now, whereas the above, exact, results are of course attributed to the precise form 
of the interpolation chosen in Eq. (14. ip for the dipole-dipole cross-section, it is clear 
that these results will approximately hold for any smooth interpolation (in between the 
shown limiting expressions) provided we impose the strong inequality Qf ^ Q2ex.p{Xy). 
Therefore, putting everything together, we arrive at (recall that y = Y — Yq) 

I pp^JXJ 1 /q2 Xy\^s 
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valid up to an overall, numerical, factor which is not under control. The effective parton 
distribution F^s of Eq. (12.91) has been generated because the saturation momenta are, 
strictly speaking, different for different types of dipoles, as explained below Eq. fl3.5l) . 
and these differences can be absorbed in the normalization of the parton distributions, as 
shown in Eq. (12. 9p . 

Apart the prefactor F^a, the above expression can be obtained from the weak-scattering 
piece in Eq. fl4.ip via the replacements r — >■ l/Qi and R — > 1/Q2- The dimensionfull factor 
I/Q2 plays the role of the 'area of the larger (target) dipole', whereas the dimensionless 
ratio Qi^'^^/Qi is recognized as the scaling variable r = Q'j.{Y)/Ql, with QsiX) the 
saturation momentum of this 'target dipole'. Accordingly, Eq. (I4.10p exhibits geometric 
scaling, as anticipated: it depends upon the resolution Ql of the 'small dipole' and the 
rapidity y only via the scaling variable r. The validity region for this behavior should be 
clear too from the previous manipulations: since the effect of the Fourier transforms is 
to select ri ~ 1/Qi and r2 ~ I/Q2 (at least, so long as the external momenta Qi and 
Q2 are well separated from each other), it is quite clear that the geometric scaling in 
the cross-section for Mueller-Navelet jets at fixed coupling holds in the same kinematical 
window as for the dipole-dipole scattering amplitude, that is. 



By the same argument, we also expect geometric scaling behavior in the 'unintegrated' 
cross-section (12.31) . for transverse momenta ki± and k2± replacing Qi and Q2 in the above 
formulae. 

5 Running coupling 

With a running coupling, the theoretical situation is less firmly under control, since 
the NLO formalism is not yet fully developed for the unitarity corrections. (The running- 
coupling version of the BK equation became available only recently [58-60].) Still, for 
the specific problem at hand, we need only some limited information about the NLO 
effects, that we believe to be reliably described by the present formalism. Indeed, to 
study geometric scaling in the Mueller-Navelet jets, we need the dipole-dipole cross- 
section in the weak scattering regime, where the BFKL approximation is expected to 
apply, and for which the NLO formalism is by now well established [11, 12] (including 
the approach towards saturation [3,7]). As mentioned in the Introduction, the use of the 
BFKL approximation is better justified in the running-coupling scenario than in the fixed- 
coupling one, since the effects of gluon-number fluctuations (which tend to invalidate this 
approximation) are drastically suppressed by the running of the coupling [52]. In fact, 
as we shall shortly discover, the main obstruction to our calculation does not come from 
our limited knowledge of the NLO perturbative formalism, but rather from a drawback 
of perturbation theory itself, ultimately associated with the running of the coupling. To 
identify this difficulty, we start by assuming perturbation theory to apply. 

As before, we choose hard cut-off momenta, Qi,Q2 ^ Aqcd, with moreover Qi ^ Q2. 
In the fixed-coupling case, this condition was enough to ensure that the relevant dipole 




(4.11) 
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sizes are sufficiently small, rj ~ 1/Qj <^ I/Aqcd, for perturbation theory to apply. With 
a running coupling, this strong correlation between Qi and rj is lost, as we shall see, but 
for the time being let us simply assume that the (effective) dipoles are perturbatively 
small. The typical situation is such that one dipole is much larger than the other, r, 
and we shall assume, once again, that the larger dipole (the 'target') evolves towards 
high gluon density and blackness on transverse sizes r < R much faster than it expands 
in impact parameter space. Under these assumptions, the dipole-dipole cross-section is 
again given by Eq. 03.21) . but with the scattering amplitude T^d{r, R,Y) now computed 
for a running coupling. This calculation has been described somewhere else [3,7,20], and 
here we present only the relevant results. 

(i) The saturation momentum is now estimated as 



Qf(i?,F) =A4cDexp 



2c(r - Yo) + pI 



P« = 111^^72 — (5-1) 



where the QCD scale Aqcd has been introduced via the running of the coupling, for which 
we used the one-loop result asiQ"^) = bo / \ia{Q^ / Aqq^)) ■ We have denoted c = feo^cAs/vr, 
with the saturation exponent ~ 4.88 (the same as in Eq. ( 13. 5p ). As before, Yq denotes 
the rapidity evolution necessary to build a saturation scale Qs equal to 1/R. 

More precisely, Eq. (15.11) has been obtained by interpolating between the asymptotic 
behavior at large Y, where the calculation is better under control, and the expected 
behavior at low Y, where one should recover the fixed-coupling result (13.61) . Indeed, for 
2cy <^ p% with y = Y ~Yq, Eq. (15. ip reduces to ~ (1/-R^) exp^Xgdy), with d evaluated 
at = 1/R^. On the other hand, for energies high enough such that 2cy ^ p^, the 
saturation momentum loses any dependence upon the target size [54]: 

Ql{R,Y) ^ Ql{Y) ^ AlcD^^ when 2cy:>p%. (5.2) 

(ii) Consider the scattering amplitude T^^^r, R, Y) for a projectile dipole with size r. 
Within the relatively wide region at 

1 « l^;:^^ ^ (2cl/ + p|)V^ (5.3) 
where the scattering is weak (Tja <^ 1), this amplitude is a universal function of the 



'scaling' variable r = r Qj{R, Y) and of Y, whose structure is quite similan I to that at 



fixed coupling, cf. Eq. (13.40 : namely, it involves the power t"^" times a function of r and 
Y which violates geometric scaling via a diffusive pattern, 
(iii) Within the more restricted window at 

1 « In^ < {2cy + piy/\ (5.4) 

the scaling violations can be neglected, and the amplitude takes the same scaling form as 
in the fixed-coupling case, that is. 
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See Eq. (3.27) in Ref. [20] for the specific function in the case of a running coupling. 
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Tdd(r,i?,F) ~ [r^Ql{R,Y)\ \ (5.5) 

Two important observations about the above results are here in order: (i) The running 
of the couphng considerably slows down the evolution, as clear from the fact that both 
the saturation momentum (15.11) and the width of the scaling region (15. 4p rise much slower 
with Y than at fixed coupling, (ii) In the high-energy regime where Eq. (15.21) applies, the 
dipole-dipole amplitude is insensitive to the target dipole size within the whole validity 
range for the BFKL approximation (cf. Eq. (15. 3p ). This second observation has dramatic 
consequences for the Mueller-Navelet process, to which we now return. 

For simplicity, we keep only the scaling piece (15. 5p in the dipole-dipole amplitude (as 
we shall later argue, our main results are independent of this approximation), and thus 
write (compare to Eq. (14. ip at fixed coupling) 

for r < rs 

fXdd(ri,r2,r) ^ 27,R' { Vt J (5.6) 

1 for r > Tg. 

where = llQs{R,Y) encodes the whole dependence upon both the rapidity and the 
target size. Below, we shall also use Tc = l/Qc(Y), cf. Eq. (15.20 . 

Let us calculate the contribution of the weak-scattering piece in Eq. (15. 6p (the first line 
there) to the Mueller-Navelet dijet cross section. (As in the fixed-coupling case, one can 
show that the respective contribution coming from the saturation piece is equal to zero.) 
We focus on the case r2 > ri, which we expect to be the most interesting one, in view of 
our condition that Qi ^ Q2- Hence, from now on, r = ri and R = r2, and we have 




■*JXJ 



F J dndr2QMQir,)Q2Ji{Q2r2)rl f^V . (5.7) 



dxidx2 

ri<rs 

It is convenient to distinguish between two regions of integration over (a) r2 > and 
(b) r2 < Tc- Given that oc exp (^—y/2cy^ is rapidly decreasing with y, it is quite clear 
(and easy to check) that the parametrically dominant contribution at large y is the one 
coming from region (a). In this region, we typically have r2 ^ Vc, so that the saturation 
scale Ts ~ Tc is independent of the size r2 of the target. (Indeed, the condition r2 ^ Vc is 
the same as p% <^ 2cy, cf. Eq. (15. 2p . applied to i? = r2.) Then Eq. (15. 7p becomes 

-T—^ ~ ^ ^ / dn Qi Ji(giri)r^ / dr2 Q2 Ji(Q2r2)r2' + . . . , (5.8) 

dXidX2 Tc" J J 

rc 

with the dots standing for the contributions coming from dipoles of size r2 < r^. We do 
our standard change of variables Ui = QiVi and U2 = Q2f^2 to obtain 

77' I ^ 



dxidx2 Qi \Qi 



(%) / du.ul^'Mui) J du2ulJi{u2) + (5.9) 



Q2/QC 
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Naturally, the regime that we are interested in is Qi ^ QcO^) ^ Q2- In this regime, QciX) 
is essentially the same as the saturation scale for a target dipole with size 1/Q2- Therefore, 
the prefactor appearing outside the integrations in Eq. fl5.9l) has the right structure to 
exhibit geometric scaling. There are, of course, additional functional dependencies in the 
limits of the remaining integrations, but at a first sight it seems that these dependencies 
are rather weak and therefore negligible: Since Qi ^ Qc, we can extend the upper limit of 
the Ml integration to oo\ then, by also making use of (12.71) . one sees that this integration 
yields a positive number of 0{1). Similarly, since Q2 -C Qc-, we can extend the lower limit 
of the U2 integration to 0. But the problem that we are facing then is that, according to 
Eq. (12. 7p . the result of the ensuing integration over U2 is exactly zero. 

Thus, by making approximations aiming at preserving the dominant contributions to 
the Mueller-Navelet cross-section at running coupling, we have found a result which is 
identically zero. Of course, a non-zero result could be instead obtained by keeping the 
formerly discarded (since formally subleading) contributions. But would that result be 
correct indeed? We do not believe so since, first, that result would be generated by phys- 
ically implausible corners of the phase-space and, second, it would be strongly sensitive 
to the fine details of our approximations — it could even oscillate between positive and 
negative values, an unacceptable feature for a cross-section. 

This invites us to critically reexamine the above calculation, in order to better under- 
stand why we obtained this vanishing result. At a mathematical level, this is related to the 
oscillatory behavior of the Bessel functions. The relevant integral, that is /q°° du2 u\ Ji(u2), 
is quite peculiar: for large values of the variable ^2, the oscillations of Ji(m2) are strongly 
amplified by the factor u\. Hence, if the overall result turns out to be zero, it is because 
of exact cancelations between large contributions with opposite signs. Then, clearly, the 
result of this integration is controlled by the behavior of its integrand at large ^2 ^ 1. If 
one sharply cuts off the integral at some value m™^^ ^ 1, then the result is an oscillating 
function of uf^^, which can vary from large positive values to large negative ones. 

We see that, even though we have tried to set up a perturbative calculation, the final, 
vanishing, result that we have obtained is in fact controlled by large dipoles — in fact, 
arbitrarily large — , for which the perturbative approach is not justified anymore. At this 
point, one may wonder about the real significance of this result: does it signal a true failure 
of perturbation theory, or is this merely an artifact of our specific approximations? The 
following argument, based on a comparison with the situation at fixed coupling, suggests 
that the first answer should be the correct one. 

The corresponding integral in the fixed-coupling case, namely /o°° d-U2 M2~^^°-^i(w2) (see 
Eq. (14.31) ). was convergent and dominated by M2 ~ 1 (i.e., r2 ~ I/Q2) because, in the 
dipole-dipole cross-section (13.21) . the rapid growth oc r| of the target area was partially 
compensated by the decay oc (l/r2)^'*''' of the scattering amplitude T^d at very large r2 
(cf. Eqs. (I3.4l) - (l3.5p with r ri and R ^ r2). In that case, the dipole-dipole amplitude 
at a given impact parameter can be made arbitrarily small be increasing the overall size 
of the target. This is a very peculiar feature of the leading-order formalism, ultimately 
related to its conformal invariance: the target size is the only dimensionfull parameter in 
the problem, so the gluon density in the target, as measured by the (local) saturation 
momentum (13.51) . must be proportional to an appropriate power of l/r2. 

The situation changes at NLO, where the running of the coupling introduces an addi- 
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tional mass scale in the problem, the 'soft' scale Aqcd- Then, for sufficiently high energy, 
the local saturation momentum becomes insensitive to the overall target size, as manifest 
on Eqs. fl5.1l) - fl5.2l) . This result is quite natural: the local gluon distribution is determined 
by the physics on the distance scale l/Qs{Y), which decreases with increasing Y , and is in 
any case much smaller than the target size r2- Similarly, the amplitude Tdd(?"i, T2, Y) for a 
small dipole (ri < l/Qs ^ ^2) and for large enough Y is insensitive to r2, as emphasized 
after Eq. flS.Sp . This property holds within the validity range (15.31) of the BFKL approxi- 
mation, and not only within the narrower window (15.41) for geometric scaling. Accordingly, 
the above conclusion about Eq. (15. 9p is more general than the geometric-scaling ansatz 
in Eq. (15.61) : within the whole range in which the BFKL approximation is expected to be 
valid, the (perturbative) dipole-dipole cross-section for two dipoles with very disparate 
sizes is expected to grow like the area ~ r| of the larger dipole. The growth is so fast that 
the cross-section for Mueller-Navelet jets is ineluctably dominated by the largest possible 
'target' dipoles, whose treatment goes beyond the scope of perturbation theory. 

Of course, in QCD dipoles cannot become arbitrarily large, because of confinement. 
In what follows we propose a heuristic modification of the previous calculation which 
limits the dipole sizes to a value ~ 1/Aqcd and thus yields a finite result for the Mueller- 
Navelet dijet cross-section. Clearly, the precise value of this result will be sensitive to 
our specific prescription for introducing confinement, and we shall try to motivate this 
prescription on physical grounds. But before we proceed, it is important to emphasize that 
this prescription will affect only the M2-integration in Eq. (15. 9p . but not also the functional 
dependencies upon Qi and F, as encoded in the prefactor there. In other terms, whatever 
prescription we choose to eliminate the large target dipoles, this will not change the 
geometric scaling behavior of the cross-section, as determined by the prefactor. 

To motivate our prescription for introducing confinement, let us recall that the 'dipoles' 
under consideration are effective dipoles, built with one gluon at x in the amplitude and 
another gluon at y in the complex conjugate amplitude, and such that R = \x — y\. 
Hence, the maximal possible value for R is the same as the maximal dispersion between 
the positions in impact parameter space at which a gluon can be produced in the proton 
wavefunction. This distance is of the order of the proton size ~ 1/Aqcd- Moreover, larger 
impact parameters ^ 1 / Aqcd lie in the tail of the proton wavefunction, where the gluon 
distribution must decay exponentially, so as it happens for any quantum-mechanical sys- 
tem with a mass gap. Thus we conjecture that the probability to produce an effective 
dipole with large size R > 1/Aqcd should fall exponentially, according to exp(— /^Aqcd)- 
We can implement this prescription in our calculation via the following replacement for 
the area factor in Eq. (15.60 (with A = Aqcd from now on): 

^ i?2exp(-i?A). (5.10) 

Then it is straightforward to see that the last integration in (15. 9p becomes 

7 1 2 / ^ / . ^ 3A 1 IQ^ 3A , , 

I d«2«2^iMexp(-.2A/g2) - g; (i + A2/g2)5/. -8g4 - g^' (^-H) 

where the approximate equality holds when A./Q2 ^ QtlQcO^) (with K/Q2 ^ 1, though). 
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a situation which is eventually reached with increasing energy at fixed Q2- Under these 
assumptions, the high-energy behavior of the Mueller-Navelet cross-section reads 



a formula which should be valid for sufficiently large Y and for Qi ^ Qc ^ Q2 ^ A. The 
overall normalization factor A/Qf in this equation depends, of course, upon our specific 
model for introducing confinement, but the scaling behavior w.r.t. r = Ql{Y)/Ql does 
not. We therefore consider this scaling behavior as a robust prediction of our analysis. 
More precisely, this behavior should hold within the window 

1 « In^^ < (2cy)l/^ (5.13) 

which is obtained after replacing r = ri 1/Qi and Qs Qc in Eq- (15.41) . (Such a 
replacement is legitimate, since the integration over ri in Eq. (15. 8p is indeed dominated 
by ri ~ l/Qi.) 



6 Conclusion and perspectives 

Our main result in this paper is that, under specific kinematical conditions — namely, 
for a sufficiently large rapidity gap Y and for a sufficiently pronounced asymmetry between 
the transverse momenta of the two jets — the partonic core of the cross-section for 
Mueller-Navelet jets should exhibit geometric scaling. 

In the leading-order formalism, where the coupling is fixed, this result is a rather 
straightforward consequence of the factorization (12.31) for the dijet cross-section together 
with known results about the dipole-dipole scattering within the framework of the BK (or 
Balitsky-JIMWLK) equation. Although the factorization (12. 3p has not been established 
here in full rigor, this should not affect the generality of our analysis, which employed 
Eq. (12.31) only in the weak scattering regime where the fc^-factorization is firmly estab- 
lished (at LO). However, as explained in the Introduction, this fixed-coupling analysis is 
a bit academic since, first, within a complete LO calculation its conclusions would be af- 
fected by particle-number fiuctuations (at least, for sufficiently high energy) and, second, 
in real QCD the coupling is running anyway. 

With a running coupling, on the other hand, our analysis lacks rigor at several points 
— it neglects other NLO corrections except for the running of the coupling and, especially, 
it turns out to transcend the framework of perturbation theory — and hence should be 
viewed as merely exploratory. Yet, we believe that our main conclusion (concerning the 
emergence of geometric scaling) is rather robust even in that context, because it is mainly 
based on the analysis of the harder jet, for which perturbation theory appears to be 
reliable. As explained in Sect. 5, the failure of perturbation theory refers merely to the 
softer jet, and, more precisely, to the connection between the transverse momentum Q2 of 
that jet and the size r2 of the associated, effective, dipole: even when Q2 is relatively hard, 
Q2 ^ -^QCD) the cross-section for dijet production is still dominated by very large 'target' 
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dipoles, with r2 ~ I/Aqcd, because the perturbative dipole-dipole cross-section grows 
very fast with r2- This growth is faster with a running couphng since the corresponding 
saturation momentum is independent of the target size. 

This failure of perturbation theory is perhaps a bit surprising, as this is not the usual 
failure associated with the BFKL 'infrared diffusion' in the presence of a running coupling: 
as expected, gluon saturation eliminates the IR diffusion and sets the argument of the 
coupling to a relatively hard scale ~ QsiY), so that the Landau pole in the coupling is 
not an issue any more. In spite of that, an infrared problem remains, as alluded to above, 
and its identification can be viewed as our second main result. 

Since based on asymptotic expansions, our results can be trusted, strictly speaking, 
only for sufficiently high energies, so their applicability to LHC may be questionable. 
This being said, and in view of the rather successful phenomenology at HERA, it would 
be nevertheless interesting to look for traces of this geometric scaling behavior in the 
forthcoming data at LHC. For instance, while keeping fixed the kinematics {Q2,V2) of the 
softest jet, one could vary the the transverse momentum cutoff Qi (with Qi > Q2) and 
the pseudo-rapidity 771 of the hardest jet, in such a way to preserve a constant value for 
the respective longitudinal momentum fraction xi ~ {Qi/ \/s)qP^ . In this way, Qi and the 
rapidity gap Y would be simultaneously changing, and then one could check whether the 
ensuing variation in the measured dijet cross-section follows indeed a geometric scaling 
pattern, cf. Eq. fl5.12p . at least approximately. Of course, such a scaling should be partially 
violated by the Qi-dependence of the parton distributions within Fofr, but this dependence 
should be rather weak and, in any case, controllable within perturbation theory. 
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